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ABSTRACT 

We observed the edge-on, planet-bearing disk of /3 Pictoris using T-ReCS at Gemini to clarify and 
extend previous observations and conclusions about this unique system. Our spectroscopy and spectral 
modeling of the 10-/xm silicate feature constrain the spatial distributions of three representative dust 
components (0.1-/xm/2.0-/mi glassy olivine and crystalline forsterite) across the inner 20-AU of the 
disk. We confirm that the 2.0-/im glassy olivine is strongly peaked in the disk center and that the 
0.1-/im glassy olivine does not show this concentration, but rather is double peaked, with the peaks 
on either side of the star. However, we do not see the strong difference in brightness between those 
two peaks reported in a previous study. Although the spatial distribution of the 0.1- /jm dust is 
consistent with the scenario of a dust-replenishing planctesimal belt embedded in the disk, we note an 
alternative interpretation that can explain the spatial distributions of the 0.1-^m and 2.0-/im grains 
simultaneously and does not require the planetesimal belt. In addition to the spectroscopy, we also 
obtained a new 11.7 /xm image of the P Pic disk. By comparing this image with that acquired in 
2003, we confirm the existence and overall shape of the dusty clump at 52 AU in the SW disk. We 
speculate that the clump's projected spatial displacement of ^2.0 AU, a 3.6-ct result, between two 
epochs separated by seven years is due to the Keplerian motion of the clump at an orbital radius of 
54.31?;!] AU. 

Subject headings: circumstellar matter — planetary systems — stars: individual(/3 Pictoris) 



1. INTRODUCTION 

The infrared emission of a debris disk arises from dust 
particles generated by collisional cascades initiated by 
large bodies like planetesimals and comets. For a small 
number of relatively nearby debris disks, it is becoming 
possible with 8-10-m-class telescopes to explore directly 
the relationship among these particles, their parent bod- 
ies, and the broader planetary cohort within a system. A 
particularly spectacular and extensively-studied example 
is P Pictoris, one of the archetypal debris-disk stars with 
which IRAS originally establish ed this important class. 

P Pic is a young (12 Myr, iZuckerman et all ,2001), 
A5 V star at a distance of 19.28 pc (jCrifo et al.lfl997li . 
The edge-on disk surrounding it shows a remarkable 
two-component structure: a main disk and a secondary 
disk that is inclined by ~5° from the main disk and 
gives rise to the inner disk warp seen in early, low- 
resolution observations ( Mouillc t et al.lU99"7t iHeap et al.l 
2000; iGolimowski et all |2006(). Although the' main disk is 
more extended and brighter in the near-IR, the secondary 
disk contributes most of the mid -IR flux and exhibits a 
conspicuous NE/S W asymmetry (jLagage fc Pantinll 19941: 
Pantin et al.lll997l ) and clumpy structures (|Telesco et al.1 
20051 here after T5). Mid-IR observations also reveal a 
central zone (r<40 AU) where dust is highly depleted 
(|Pantin et alj |1997). 10-^m-band spectroscopy of P Pic 
shows the silicate feature similar to that of comet Halley, 
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which suggests a cometary origin of the dust in the in- 
ner di sk of P Pic (iTelesco fe Knackdll99"T[ iKnacke et aLl 
II991 lAitken et alJH993D~ This conclusion is supported 
by the modeling by iLi fc Greenberel (119981) and obser - 
vations at other mid-IR wavelengths (jChen et al.ll2007T) . 
However, collisional processes among km-sized planetes- 
imals may also account for the observed dus t properties 
(jThebault et all 120031 IGrigorieva et al.ll2007t ). 

Given the wealth of observational constraints, con- 
siderable effort has been expended to produce a co- 
herent picture of the P Pic s ystem and properties 
of it s dust (e.g.. IPantin et al. 19971: ILi fc Greenbersj 
1991 iThebault et all 120031: IGrigorieva et all 120071: 
Czechowski fc Mannll2007t lAhmic et alll2009l ). It is gen- 
erally agreed that dust grains are being produced con- 
tinuously by comets orbiting close to the star, collisional 
cascades, and/or the collisional break-up of planetesi- 
mals, with the particle evolution being strongly affected 
by the stellar radiation and the gravity of planets in the 
disk. 

Some features of the disk of p Pic, such as the in- 
clined secondary disk and its asymmetric morphology, 
suggest the presence of a planet; in fact, one was dis- 
covered recently by direct i maging (|Lagrange et al.|[2009l 
l20l0llBonnefov et alJ2011l) . The most recent and precise 
measurements of the position angles (PA) of the main 
and secondary disks and P Pic b imply that the planet 
is located in the se condary disk (|Lagrange et all 120121 : 
IChauvin et alll2012D . thus providing a promising expla- 
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nation for it s origin. However, at present, other possi- 
bilities (e.g.. lGrigorieva et al.ll2007h cannot be excluded, 
and new clues may illuminate the relationship between 
the planetary system of ft Pic and the observed disk mor- 
phology. 

10-Mm-band spe ctroscopic observations by 
IQkamoto et all (l200l hereafter 04) using COMICS 
at the Subaru Telescope reveal what may be several 
planctesimal belts in the ft Pic disk, the most prominent 
of which is located at 6.4 AU from the star and possibly 
tilted with respect to the secondary di sk. T5 observed ft 
Pic using T-ReCS (jTelesco et al.lll998f > at Gemini South 
in 2003, and their multi- wavelength images showed a 
bright clump in the disk at about ~52 AU southwest 
of the star and composed of particles that may differ 
in size and/or composition from dust elsewhere in the 
secondary disk. This bright clump of thermally emitting 
dust may mark the location of collisionally grinding, 
resonantly trapped planetesimals or, as T5 emphasized, 
the cataclysmic break-up of a planetesimal. 

To clarify and extend some of these previous observa- 
tions and conclusions about the ft Pic disk, we present 
in this paper new 8-13 /tm spectra and 11.6 /im (Si-5) 
images made with sub-arcsec angular resolution across 
the 40-AU-diameter and the 160-AU-diameter regions in 
the spectroscopy and imaging, respectively. While the 
spectra by 04 were obtained at high airmass (>3.0) and 
suffered from strong atmospheric ozone absorption near 
9.8 /im, our new data, made recently using T-ReCS at 
Gemini South, were obtained under more favorable con- 
ditions and so deal with some ambiguities in the 04 ob- 
servations. Using T-ReCS, we also repeated a subset of 
observations by T5. We use these two-epoch observa- 
tions, both made with T-ReCS and spanning a period 
of ~7 years, to constrain further the properties of the 
bright clump in the disk of ft Pic. 

The details of our observations and data reduction are 
presented in the following section. The analysis and mod- 
eling of the data are covered in the Section 3. Finally, 
Sections 4 and 5 are devoted to discussions and conclu- 
sions, respectively. 

2. OBSERVATION AND DATA REDUCTION 
2.1. Spectroscopy 

Spectra of ft Pic were obtained during 2010 December 
16-18 (UT) with the low-resolution-10-/im mode (nomi- 
nal resolution R = 116 at 10 /on) of T-ReCS. A 0.35" 
slit was oriented along the disk at position angle 32° , the 
disk inclination measured by T5 using N- and Q-band 
images. Note that this value is very close to, but not ex- 
actly, the orientation o f the secondary disk m easured in 
the near-IR (32.7-33.8°. [Lagrange et"atll2012| ). The pixel 
scale was 0.09". Eight integration sets of ft Pic were ac- 
quired to give a total on-source time of 2,080 s, which 
was about 60% of the time awarded for the program. 
The science observations were interlaced with those of 
the Cohen standard HD 39523 (Kl III, 5.1° away from ft 
Pic) in order to monitor the point spread function (PSF) 
and the flux/telluric correction. All data were taken in 
the standard chop-nod mode with a chop frequency of 
2.3 Hz and a chop throw of 15" perpendicular to the slit. 

These data were reduced with custom IDL routines. 
All frames in the raw T-ReCS FITS files were checked 



Table 1 

Summary of T-ReCS imaging of /3 Pic 



2003 



2010 



Program ID 
UT Time 

Integration Time (s) 
Filter 

Pixel Scale (arcsec) 
Mean Air Mass 
FWHM a (arcsec) 
Instrument PA (degree) 
Chop PA (degree) 



GS-2003B-Q-14 GS-2010B-Q-50 
2003 Dec 30 2010 Dec 16 

2 x 450 3 x 300 

Si-5 (11.7 fim) 
0.09 

1.196 1.074 
0.36 0.35 
340 
123 122 



a Measured in the PSF references. 

visually and ^25% of them with low quality (e.g., bad 
seeing, strong background) were discarded before fur- 
ther reduction. Other artifacts, like non-flat background, 
channel offsets, jitter, and bad pixels, were corrected in 
each FITS file before combining the data. 

The wavelength calibration was done using the telluric 
lines identified in the raw spectra. The flux and telluric 
corrections were performed using the Cohen m odel tem- 
plate spectrum of HD 39523 (jCohen et al.lll999D . We also 
calibrated the absolute flux at 11.7 /im using our imag- 
ing data, and these two methods agreed very well (within 
the typical mid-IR photometric accuracy of 10%). 

One spectrum was acquired for each spatial pixel along 
the slit (corresponding to a 1.74-AU-wide region at the 
distance of ft Pic) within 20 AU to the star. As a 
double check of the photometric calibration accuracy, 
we integrated all reduced spectra within an aperture of 
3.7" centering on the star to mimic a single, low-spatial- 
resolution spectrum of ft Pic, and the r esult matched 
perfectly the data by iKnacke et al.l (|1993l) . 

In order to separate the disk excess from the total 
flux (star plus disk) at each wavelength, the spectral en- 
ergy distribution (SED) of the stellar photosphere was 
assumed to be that of an 8,200 K blackbody with a flux 
de nsity of 2.85 Jy at 8. /im. This value was estimated 
by IKnacke et al.l Jl993) by matching the photospheric 
continuum to the near-IR photometric measurements of 
ft Pic. The uncertainty of their fit was not stated, but 
as we show in the following section, the uncertainty in 
the stellar flux should not affect the conclusions of our 
spectral analyses. The PSF profile at each wavelength 
was assumed to be Gaussian, the FWHM of which was 
0.5", measured from the spectrum of HD 39523. 

2.2. Imaging 

The two-epoch T-ReCS imaging data were obtained on 
2003 December 30 (T5) and 2010 December 16 through 
the same Si-5 (11.7 /im ) filter. Key information about 
these two observations is summarized in Table [U Final 
images were co-added using frame registration to achieve 
the highest spatial resolution that the raw data could 
deliver. 

To minimize the introduction of any spurious fea- 
tures due to row- or column-correlated pattern structure 
known to be associated with the detector array, T-ReCS 
was rotated so that the disks major axis was tilted in the 
raw images. During the reduction, we rotated the data 
again to display the disk horizontally in the reduced im- 
ages. 
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Because our major science focus with the images is to 
examine the morphology, rather than the total flux, of 
the bright clump in the disk, we did not use the Cohen 
standard for the calibration. Instead, both images were 
normalized by the peak intensity of /3 Pic measured in 
each image through the Lorentzian profile fit. 

3. RESULTS AND ANALYSES 

3.1. Spectroscopic decomposition 

The reduced N-band spectra (Figure [T]) show clearly 
the broad 10-/im emission with a shape that depends on 
position across the inner disk of (3 Pic. The shape of the 
10-/im silicate feature, often quantified by the flux ratio 
of -F11.3/-F9.8, has long been considered a measure of the 
degree of dust processing (dust growth and crystalliza- 
tion, e.g.. iBouwman et all 120011 : Ivan Boekel et all [20051: 
iManoi et al.ll201lD . Our spatially resolved spectra can be 
used to constrain the spatial distributions of amorphous 
and crystalline silicates across the disk of f3 Pic, thus 
providing information on the processing and evolution of 
the dust. However, due to the strong telluric ozone fea- 
ture around 9.8 /im, obtaining a reliable estimate of the 
-F11.3/-F9.8 ratio is a challenge. In order to mitigate this 
pr oblem, we followed a n alternative approach developed 
bv IHonda et al.1 (|2003h and 04, and fit the entire 10-/xm 
silicate feature with a spectral model composed of four 
emission components: 0.1-/xm and 2.0-/im glassy olivine, 
sub-micron-sized crystalline forsterite, and a power-law 
continuum. The combination of the first two compounds 
was considered representative of a large size range (0.01- 
5 /im) of amorphous silicates th at are thought to dom - 
inate the 10-/zm silicate feature (Bouwman et^l200l . 
The power-law continuum mimics all relatively feature- 
less dust emission in the 10-/im window. The fitting for- 
mula is 



= a Q (X/10fim) n +Y^ fliK'i(A)B (A, 7*) /B (10/im, T,) . 



Here 
W cm~ 2 /im -1 



is the observed brightness (in 
arcsec -2 ) at wavelength A (in /im), 
a is the normalization factor, and B is the Planck 
function. The parameter n defines the slope of the 
power-law continuum. The subscript i indicates four 
spectral components: continuum {i = 0), 0.1-/im 
glassy olivine (i = 1), 2.0-/im glassy olivine (i = 2), 
and crystalline forsterite (i — 3). The parameters k\ 
are the normalized mass absorption coefficients, and 
therefore are dimensionle ss. The normalized factor for 
each species is given by iHonda et al.l J2003). For the 
purpose of comparison, we chose to use the same data 
set of tabulated values of k as used in 04. For the 0.1- 
and 2.0 -/im glassy olivine, k is calculated using Mie 
theory (jDorschner et al.l 119951: IHonda et all l2003fi: for 
the cr ystalline forsterite, k is measured by iKoike et al.l 
(2003). T is the dust temperature (in K). According 
to such a fitting formula, any remnant photospheric 
emission (due to inaccurate stellar subtraction) is not 
distinguishable from the underlying continuum, i.e., the 
first term in Equation 1, so it should not affect the 
fitting results of the other three dust components. 
As proposed by some studies, the mid-IR spectrum 
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Figure 1. Some example spectra of /3 Pic (thin solid lines), 
overlapped with the best-fit spectral models (thick lines). Spec- 
tra are disk emission only, i.e., with stellar photospheric emission 
subtracted. Model spectra consist of four components: 2.0-^tm 
glassy olivine (dashed), 0.1-/im glassy olivine (dotted), crystalline 
forsterite (dash-dotted), and power-law continuum (not shown). 
Error bars for each fitted component are shown in Figure [2] 

of /3 Pic can be better reproduced by assuming core- 
mantle silicate particles (i.e., a cometary dust model), 
rather than the popular Draine and Lee interstellar sil- 
icates (e.g iDraine fc Led 119841: iLi fc Greenbergl 119981 : 
IChen et al.ll2007D . However, we find that the latter fit our 
data fairly well, so we do not consider the core-mantled 
silicates in our work. 

For amorphous olivine grains (i.e., i = 1, 2), the tem- 
perature is derived by balancing the radiative absorption 
and re-emission (assuming spherical grains): 



«i (A) ( — 



B x {T,)dX = / An t {X)B x (T t {r))dX 
Jo 

(2) 

where i?* and T* are the stellar radius (6.8 x 10 3 AU) 
and effective temperature (8,200 K) of j3 Pic. The param- 
eter r is the disk radius, in AU, at the center of each spa- 
tial pixel. The temperature of the crystalline forsterite 
grains is assumed to be 1,000 K, i.e. approximately their 
annealing temperature. We note that, however, the fit- 
ted brightness distribution of forsterite is not sensitive 
to this assumption. Assuming a higher (e.g., 1,200 K, 
as in 04), lower, or varying temperature (i.e., similar to 
the small and big glassy olivine) affects the final results 
little. 

For each spectrum associated with a certain spatial 
pixel, there are five free parameters, a, (i — 0,1,2,3) 
and n, to fit the spectral model. The best fit was ob- 
tained by minimizing the reduced chi-square (xt) over a 
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a 0.1 um glassy olivine x 2.0 
■ 2.0 um glassy olivine 
x Crystalline forsterite x 2.5 
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Figure 2. Emission distributions of small (open triangles), big 
(filled squares), and crystalline grains (crosses) in the fj Pic disk. 
Note that the profiles are multiplied by different factors for the pur- 
pose of illustration. The approximate spatial resolution (FWHM) 
of T-ReCS at 10 fim is indicated in the top-right corner of the plot. 



reasonably large parameter space. The goodness of fit 
was also checked by eye, and we determined that all fits 
for which the xt deviated from the minimum xt by more 
than 10% were visually distinguishable from the best-fit 
model, and were considered "bad fits" . In contrast, those 
for which the deviations were less than 10% were consid- 
ered "good fits" . For each free parameter, we then used 
all the values given by those "good fits", and defined 
the standard deviation of that ensemble of values as our 
stated 1-ct uncertainty for that parameter. 

Figure [TJ shows some examples of the model fit results. 
The data quality allows us to extract the 10-/xm silicate 
feature out to about 20 AU (~T") from the central star. 
Note that the ozone-dominated wavelengths (9.4-9.8 /im) 
were not used in the evaluation of xt- Although the peak 
around 11.3 /im does not appear to be particularly strong 
in the spectra, we found that excluding the crystalline 
forsterite component increases the xl by at least 20% 
for the central pixels where the signal-to-noise ratios are 
high, and therefore we decided to keep this component 
in the spectral decomposition and analysis. 

The spatial profiles of ai, a-i-, and 03 are plotted in Fig- 
ure [2] For a direct comparison with the result presented 
by 04, the spatial profiles are binned every two pixels 
(except for the central one) to closely match the spatial 
sampling used in 04. Figure [5] shows clearly that the 
brightness distributions of three dust species are differ- 
ent from each other, as concluded by 04. The 2.0-/zm 
glassy olivine particles (squares) have a single-peak dis- 
tribution centered on the star; their brightness profile 
has a FWHM of ~7 AU, comparable to the measured 
FWHM of the spectrum of the standard star (~9 AU 
at 10 /im). In contrast, the 0.1-/zm glassy olivine par- 
ticles (triangles) show two almost symmetrically placed 
peaks, each at ~5-6 AU from the star and ~8 AU wide. 
The profile of crystalline forsterite particles (crosses) is 
weakly centrally peaked and symmetric within the un- 
certainties. We conclude that the crystalline forsterite 
particles are distributed across a region comparable in 
scale to that of the 2.0- /im glassy particles and unlike 
the double-peaked 2.0-/im glass olivine particles. 
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Figure 3. Smoothed Si-5 (11.7 M m ) images of /3 Pic, obtained 
in 2003 (upper panel) and 2010 (lower panel); NE to left, SW to 
right. The hatched circle in the lower left of each panel represents 
the spatial resolution achieved with these observations, being the 
smoothed FWHM of a point source. Both images have been rotated 
counterclockwise by 48° (2003), or 68° (2010). The vertical lines 
are at star (dotted) and at SW 52 AU (dashed). Contours are 
linearly spaced by 3-cr, starting at 3-cr. The features labeled by 
crosses are not real but ghost images of Pic. 

3.2. Imaging 

Here we compare the mid-infrared images (smoothed 
with a Gaussian kernel of 3-pixel in FWHM) of ft Pic 
obtained in 2003 and 2010, focusing on the bright clump 
at 52 AU in the SW wing of the disk. For a compre- 
hensive analysis of the 2003 images, per se, we refer the 
reader to T5. Two-epoch data are shown together in 
Figure Note that some small morphological changes 
associated with the different instrument setups are evi- 
dent. One of them is the ghost image of ft Pic (labeled 
by crosses), which rotated by ~ 20° with respect to ft 
Pic between two observations. This value is consistent 
with the change in instrument PA (Table [TJ . Neverthe- 
less, the two data sets were obtained under similar sky 
conditions, and their image qualities are comparable in 
terms of their FWHM and signal-to-noise ratios. 

The bright clump, discovered by T5, embedded in the 
SW wing is unambiguously detected and resolved in both 
images. Note that the proper motion of ft Pic is ~ 0.5" 
(i.e., 9.6 AU at the distance of 19.28 pc, based on Hippar- 
cos data. Ivan Leeuw en 2007) over seven years, a displace- 
ment that would have been obvious in the two-epoch 
images if the clump were a chance superposition of an 
unrelated background object; this observation confirms 
the conclusion of T5 that the clump is associated with 
the disk, and gives rise to the SW/NE asymmetry ob- 
served in the disk in earlier mid-IR observations (e.g., 
iPantin et aLlll997ft . 

Aided by the line drawn at SW 52 AU, we can see 
that there seems to be a small offset between the 2003 
and 2010 locations of the clump. As was done by T5 
to investigate this region more carefully, we subtract the 
emission in the fainter NE wing from the SW wing, and 
draw the contours of the residual emission (assumed to be 
mainly contributed by the bright clump) in the left panels 
of Figured] The profiles plotted in the lower- right panel 
show the flux distributions along the major axis of the 
disk. As a comparison, the best-fit Lorentzian profiles of 
ft Pic are also plotted in the upper-right panel, where the 
two-epoch data match each other almost perfectly; thus, 
there is no significant change in the PSF morphology be- 
tween the 2003 and 2010 observations. Furthermore, if 
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Figure 4. The residual emission (left panels) of the SW 52 clump, 
resulting from subtracting the NE wing from the SW wing along 
a line through the star. Contours are spaced by 3-cr, starting at 
3-cr. The profiles plotted in the lower-right panel show the varying 
flux (normalized by the peak intensity of j3 Pic) along the major 
axis of the disk. The cross drawn at right shows the smoothed 
point source FWHM (horizontal bar) and the photometric error 
(vertical bar). The alignment uncertainty between the 2003 and 
2010 profiles is estimated to be well below 1 pixel (i.e., 1.7 AU at 
the distance of Pic). The profiles drawn in the upper-right panel 
are best-fit Lorentzian PSF models of f) Pic obtained in the same 
data set. 

we divide each data set of 2003 or 2010 into two sub- 
sets differing in observing time and repeat the reduction 
procedures, the offset is still visible, which implies that 
changes in seeing are unlikely the cause of the offset. 

To quantify the spatial displacement shown in the 
lower-right panel of FigureSJ we performed a least-square 
fitting of the Lorentzian PSF model to the vicinity (—18 
AU in width) around the clump in each image. The best- 
fit profiles suggested a radial shift of 2.0 AU between the 
2003 and 2010 data. To estimate the error, we consid- 
ered the following sources of uncertainties: (a) the cen- 
tering error, i.e., how accurately the two-epoch data were 
aligned using the position of (3 Pic; (b) uncertainty in the 
relative image orientations; (c) image distortion intro- 
duced by the image rotation algorithm. For (a) and (c), 
we performed a series of tests with artificial stars, and 
determined that under a typical noise level (S/N — 1000), 
a 0.2-pixel centering accuracy is achievable; the distor- 
tion due to the image rotation algorithm is <C0.1 pixel 
and thus negligible. We also checked the commissioning 
documents of T-ReCS, and found that the amount of im- 
age distortion is linearly proportional to the field of view. 
The maximum distortion is about 1 pixel, and occurs at 
four corners of the detector array. Since the disk of j3 
Pic occupies only a small portion around the array cen- 
ter, and the relative image orientation between the 2003 
and 2010 data is merely 20°, we estimated that the error 
due to (c) is no greater than 0.25 pixel. Taking all factors 
discussed here into account, the total uncertainty, which 
is the quadratic sum of (a), (b), and (c), is 0.32 pixel 
(i.e., 0.55 AU). Based on this, we conclude that the spa- 
tial displacement of the clump between 2003 and 2010 is 
marginally significant (3.6-cr) and therefore worth some 
consideration at this time. However, additional obser- 
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Figure 5. Comparison between our work and the 04 results on 
the spatial distributions of the 0.l-fim/2.0-fim glassy olivine. Note 
that different normalization factors are used for 04 and our data 
to aid visual comparison. 

vations must confirm it over the coming years. We will 
discuss its scientific implications in the following section. 

4. DISCUSSION 
4.1. Radial distributions of dust grains 

The spatial distributions of the three representative 
species derived by our spectral modeling have both sim- 
ilarities and important differences with those presented 
by 04. As seen in Figure [5j both studies show that the 
2.0-/im glassy olivine is strongly peaked in the disk cen- 
ter, while the 0.1-/im glassy olivine does not show this 
concentration, but rather is double peaked. There is an 
apparent asymmetry in the brightness profile, and hence 
spatial distribution, of the 0.1-/im glassy olivine compo- 
nent. However, it is not as strong as that inferred by 04, 
and its significance not as firm. Also, we see a weaker 
crystalline feature: the total emission of forsterite (i.e., 
the profile shown in Figure [5] integrated along the entire 
disk) is —20% of that of the 2.0-/im glassy olivine, in 
contrast to —40% in 04. 

In the two studies (04 and ours), the slits were placed 
along different PAs (30° and 32°, respectively), but this 
is not sufficient to explain the disagreement. The slits 
used in both studies are —0.35" wide, thus not particu- 
larly narrow compared to the slit length of —1.5" (i.e., 
the diameter of the disk where the silicate feature is 
clearly detected). The offset due to the 2° difference in 
PA is much smaller than the slit width, and thus should 
have little effect on the observations. Another difference 
between the two studies is the way the dust tempera- 
ture is calculated. We found that the equation used by 
04 is only appropriate for particles larger than typical 
wavelengths of absorbed stellar radiation, an assumption 
that is not valid for 0.1-/im grains, so we switched to an- 
other method as described in Section 13.11 However, we 
could not reproduce such a strong asymmetry, even if we 
adopted the same equation as in 04 to calculate the dust 
temperature. 

We speculate that the different degrees of asymmetry 
observed for 0.1-/im grains by us and 04 can be explained 
by either a dynamical evolution of the inner disk of /3 Pic 
during the time interval between our observations and 
those of 04, or an observational effect associated with 
the positioning of the slit by 04. At a radius of 5 AU, 
the dynamical timescale (i.e., the orbital period) in the 
disk of /? Pic is comparable to the time span between 
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two observations (2003 and 2010, respectively). Thus, 
the change in the degree of asymmetry may be caused 
by the orbital movement of the parent dust-replenishing 
bodies that are not evenly distributed in the belt. How- 
ever, 04 suggest that the asymmetry may be at least 
parly attributable to a slight de-centering, or shift, of 
the slit position toward the SW from optimum. In addi- 
tion, the ozone feature between 9.4-9.8 /im can make the 
fit of the 0.1-/mi glassy olivine problematic, especially at 
higher airmasses such as those (>3) which apply to the 
04 data. Finally, the discrete spatial sampling associated 
with pixelation (i.e., the disk of j3 Pic not being exactly 
symmetrically sampled by the pixels along the slit) may 
contribute to the observed asymmetries. This last effect 
is less important for our data, because the T-ReCS pix- 
els are smaller. In the end, however, we conclude that 
there is currently no definitive way to account for the 
differences in our observations and those of 04. 

4.2. Crystalline forsterite 

As seen in Figure [TJ the crystalline emission peak 
at 11.2 /im is not strong in our spectra (cf. Figure 1 
of 04). The resultant spatial distribution of the crys- 
talline forsterite is mildly centrally peaked, being several 
times broader at FWHM than that of the 2.0-yum glassy 
olivine. This result is consistent with the scenario of 
dust crystallization through thermal annealing by stel- 
lar radiation. In the disk of /? Pic, micron-sized grains 
can reach the classical annealing temperature of 1,000 
K at ~0.6 AU from the star. At such temperatures, 
the annealing timescal e (of the order of hours to days, 
Hallcnbec k et al.l [1998) is much shorter than the colli- 
sional timescale (> 10 4 - 10 5 years, lAhmic et al.|[2009h . 
In order to form an observable concentration, the newly 
formed crystalline grains must be orbitally stable, i.e., 
with a ratio of radiation pressure to gravity < 0.5. This 
holds true for grains larger than ~1.3 /xm (jArtvmowiczl 
1997). Much smaller (e.g., O.T/zm) grains can reach 
the annealing temperature at larger distances from the 
star, but are more affected by radiation pressure and 
repelled quickly (orbital timescales) from the star. For 
intermediate-sized grains, i.e., those with /3 close to the 
blow-out limit, some degree of radial mixing following 
the thermal annealing is probable, and may account for 
the broadening in the spatial profile of forsterite rela- 
tive to the 2.0-//m glassy particles. We conclude that, 
given the similar spatial distributions and characteristic 
particle sizes, it is possible that the crystalline forsterite 
grains are the annealed descendants of the larger glassy 
olivine particles. 

Our admittedly speculative scenario clearly does not 
exclude other mechanisms that have been proposed to 
explain the formation of crystalline silicate minerals and 
their existence in the cold, out er disk. Those mech- 
anisms, such as shoc k heating (jNuth fc Johnson! 120061 : 
I Alexander et alJl2007h . aqueous alteration of dust on as- 
teroidal parent bodies (jNuth et al.l I2005D . and crystal- 
lization associated with the formation and subsequent de- 
struction of mas sive fragments in young protostellar disks 
(|Vorobvovl2011h . can occur at large distances rather than 
the inner few AUs, but their signatures may not be ob- 
servable by us. 

4.3. A planetesimal belt 



The double-peaked spatial profile of the O.T/Ltm grains 
is interpreted by 04 as tracing a planetesimal belt that 
continuously replenishes the disk with these amorphous 
silicate grains. Due to the high /3-value, sub-micron 
grains should be blown out of the disk in a very short 
times cale (less than several tens of years, lArtvmowiczl 
1988), and therefore the peaks of their spatial distri- 
bution should trace the dust-replenishment sites. Our 
observations imply that the radius of the belt proposed 
by 04 as the origin of the 0.1-/im grains is 5.6 ± 3.0 AU 
(assuming Gaussian profiles for each "lobe" ) . If the plan- 
etesimal belt is an analogue of the asteroid belt in the 
Solar System (i.e., the belt is formed by the gravitational 
perturbation induced by a giant planet, and the inner 
and outer boundary of the belt are defined by the 4: 1 and 
2:1 resonance with the planet, iFerraz-Mellol 119941 ) . the 
measured location of the belt implies that the perturbing 
planet has an orbital radius of 9-14 AU. This is in rough 
agreement with the most up-to-da te measurement of th e 
orbital radius of [3 Pic b (8-9 AU, iChauvin et "aTll2012ft . 
In addition, the orientation of the orbita l plane of Pic 
b wa s recently determined to be ~212° ((Lagrang e et al.l 
I2012T ) . i.e., in good alignment with the inner disk orien- 
tation measured in the mid-IR (T5). 

4.4. An alternative interpretation for the different 
distributions of small/big grains 

As already noted, the brightness distribution of the 
2.0-/im grains is very different from the double-peaked 
profile of the 0.1-/xm particles. 04 account for this dif- 
ference by proposing that the observed 2.0-/j,m grains, 
which they assume originate from the same planetesi- 
mal belt as the O.T/xm particles, are propelled towards 
the central star by Poynting-Robertson drag, eventually 
concentrating near the disk center. In fact, this mech- 
anism probably does not work in a disk as massive as 
/3 Pic, where the collisional timescale is t hought to be 
much shorter than that of the P-R drag (|Artvmowiczl 
119971: lWvattll2005l see also lAhmic et al.ll2009fl Rather, 
most of the grains do not travel far from their birth- 
places before mutual collisions grind them fine enough 
to be expelled by radiation pressure. Therefore, we pro- 
pose an alternative explanation for the difference in the 
brightness distributions of the two grain components, as 
follows. 

Considering that small grains can reach a significantly 
higher temperature than do larger grains at the same 
radius of the disk (as shown in Figure [6j where tem- 
perature is calculated using Equation 2), the remarkable 
spatial profile of the OA-fim glassy olivine may be ex- 
plained by a strong depletion of such grains around the 
disk center, as a consequence of thermal processing (an- 
nealing, sublimation) of the dust. That is, those par- 
ticles at the center are removed or destroyed. Accord- 
ing to Figure [51 0.1-/im glassy olivine at ~3 AU can be 
heate d to 800 K, for which the annealing tim escale is ~ 
3 yr (jFabian et al.l 120001: ITanaka et al.l 120101) . i.e., com- 
parable to the radiation pressure blow-out timescale (for 
sub-micron grains) and the orbital period at the radius 
of several AU, but significantly shorter th an the colli- 
sional timescale of sever al thousand years (|Artvmowiczl 
[19971 lAhmic et al.l[2009h . Since the annealing and subli- 
mation timescales decrease abruptly as temperature in- 
creases (|Tanaka et al J 12010). a deep depletion of 0.1-/im 
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Figure 6. Dust temperature calculated using Equation 2 for 0.1- 
)im and 2.0-fim glassy olivine, assuming a stellar luminosity of 8.7 
Lp). Annealing and su blimation temperatures are adopted from 
IHenning Meeusll200g| . 

amorphous grains can be expected not far inside this 
radius. To reach the same temperature, 2.0-/zm grains 
need to be as close as 0.6 AU to the star, and the cen- 
tral depletion zone would be too narrow to be detected 
by our observations. Thus, in this scenario, both dust 
species originate from parent dust populations that are 
centrally peaked. Their steady state distributions in the 
disk are modified by thermal annealing and sublimation, 
and coarsely sampled by our observations. 

To see whether a central clearing zone of diameter 6 AU 
for the small amorphous grains is even resolvable by T- 
ReCS, we construct a toy model of the disk, in which the 
spatial distribution of 0.1-/im grains is Gaussian shaped 
but modified so that is has a central sharp-edged hole of 
radius ri n . This corresponds to a particle surface den- 
sity peaked near rj„ and then decreasing gradually out- 
wards. We also assume a uniform dust temperature and 
optically thin disk, so that the emission brightness is di- 
rectly proportional to the amount of dust encountered by 
the line-of-sight. Finally, the intrinsic brightness profile 
is convolved with the PSF to mimic the real observa- 
tions. In Figure the results are plotted against the 
observed spatial profile of the 0.1-/Ltm glassy olivine. In 
all simulations, the FWHM of the PSF is set to be 7 
AU, i.e., roughly the width of the brightness profile of 
2.0-/im grains. We find that resolving a central clearing 
zone with ri n of 3 AU would be a challenge for T-RcCS. 
Nevertheless, the hole would be clearly seen when ri n is 
around 4-5 AU, close to the value of ri n suggested by 
the dust annealing temperature as discussed above. We 
also note that our model is unable to reproduce the great 
depth of the central dip, which is very noticeable in the 
real data but not so in the simulations. This may indicate 
that the seeing was significantly better during the obser- 
vations of j3 Pic than when the PSF was observed, or it 
may reflect additional uncertainties in the photospheric 
subtraction or the spectral fitting. 

4.5. The dust clump at SW 52 AU 

Using multi-wavelength mid-IR images, T5 explored 
the properties of the clump at SW 52 AU. The SED of the 
clump implies a characteristic particle temperature there 
of ^190 K, noticeably higher than the value of about 140 
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Figure 7. Brightness profiles of the 0.1- fim glassy olivine, derived 
from real observations (solid line) or simulations (dotted lines), in 
which the r^ n (radius of the central clearing zone) are assumed to 
be (from bottom to top) 3, 4, 5, 6, and 7 AU, respectively. In all 
simulations, the spatial resolution is assume to be 6 AU. 

K for the region on the opposite side of the disk, at NE 
52 AU. Therefore, T5 infer that a population of sub- 
micron-sized (0.1-0.2 fj,xn) particles exists in the clump, 
with their sizes and/or compositions different from the 
dust elsewhere in the disk. They conclude that the catas- 
trophic collisional break-up of a planetesimal (~100 km 
or larger in diameter) is a plausible explanation for this 
localized inhomogeneity in the dust population. In this 
picture, fragment post-collision velocities of ~l-2 km s _1 
and radiation-pressure-driven dispersal account for the 
vertical and radial extent of the clump. Several stud- 
ies (e.g.. iDermott et all 120021 : iKenvon fc Bromleyl l2004t 
IGrigorieva et al.l 12007ft indicate that such catastrophic 
collisions should indeed be observable. 

The two images taken seven years apart may shed ad- 
ditional light on the clump properties. Figure [8] shows 
the projected displacement over seven years expected for 
an object orbiting /3 Pic in an edge-on, circular Keple- 
rian orbit. If we assume that the orbital movement of 
the clump is traced by the locations of the peaks in the 
emission residuals shown in Figure [4j then the measured 
projected shift of 2.0 AU between 2003 and 2010 is con- 
sistent with a true orbital radius of 54.31^2 AU. Since 
the SW disk is rotating toward us, as su ggested by high 
resolu tion Na I D spectroscopy of j3 Pic (jOlofsson et all 
120011) . the observed shift also implies that clump is mov- 
ing towards us, and approaching its greatest southwest- 
ern elongation. This tentative conclusion can be tested 
over the coming years, and, if confirmed, would consti- 
tute unique insight into the dynamics of a debris disk; to 
date, e Eri is the only other debris disk fo r which orbital 
prope r motions may have been measured (|Greaves et al.l 
120051 ). 

Also note that some features of the dust clump, like 
a high proportion of sub-micron grains, physical exten- 
sions in both vertical and radial directions, and the or- 
bital movement, are in general consistent with the con- 
sequence of a collisional avalanche, i.e., a chain reaction 
initiated by a break-up of a comet or a planetesimal-belt 
object, producing a substantial amount of small grains 
propagating outwards in th e debris disk (|Artvmowiczl 
Il996t IGrigorieva et alj |2007). Although in this scenario, 
the current morphology of the clump indicates that we 
are witnessing a rare event only several tens to several 
hundred year after the catastrophic break-up, such a pos- 
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Figure 8. The projected displacement expected for an object 
orbiting f) Pic in a circular, edge-on Keplerian orbit over 7 years, 
assuming a stellar mass of 1.75 Mq. The vertical error bar indicates 
the observed radial displacement and its uncertainty (2.0±0.6 AU), 
corresponding to the most probable orbital radius of 54.3^j'j AU. 

sibility cannot be ruled out. Current data are not sen- 
sitive enough to allow us to explore the grain properties 
in this region as we have done for the innermost part of 
the disk. Deeper spectroscopic observations in the future 
will be crucial to understand the nature of this striking 
feature in the disk of j3 Pic. 

Finally, we note that, although the clump is clearly 
seen in the mid-IR, no counter part has been found at 
optical or near-IR wa velengths (Go limowski et aD 120061 : 
iBoccaletti et al J 12009( 1. This is not surprising, since the 
scattered light and mid-IR emission may trace different 
regions of the disk. Nevertheless, this issue has not been 
studied in detail, and important hints of the dusty envi- 
ronment and formation mechanism of the clump can be 
revealed by a model that is able to reproduce the multi- 
wavelength observations simultaneously. 

5. SUMMARY 

We observed the edge-on, planet-bearing disk of /3 Pic 
using T-ReCS at Gemini to clarify and extend previous 
observations and conclusions about this unique system. 
We used spectroscopy and spectral modeling of the 10- 
/im silicate feature to constrain the spatial distributions 
of three representative dust components (0.1-/im/2.0-/im 
glassy olivine and crystalline forsterite) across the inner 
20- AU of the disk. We also obtained N-band images of 
the /3-Pic disk. We conclude the following: 

1. We confirm the 04 results that the 2.0-fxm glassy 
olivine is strongly peaked in the disk center and that, in 
contrast, the 0.1-/mi glassy olivine does not show this 
concentration, but rather is double peaked, with the 
peaks on either side of the star. However, we did not 
see the strong brightness asymmetry between those two 
peaks of the 0.1-/xm glassy olivine reported by 04. 

2. Compared to the observations by 04, the distribu- 
tion of the crystalline silicates observed by us is much 
weaker and not so strongly centrally peaked. 

3. Although the spatial distribution of the 0.1-/zm 
glassy olivine is consistent with the general scenario of 
a dust-replenishing planetesimal belt embedded in the 
disk, as proposed by 04, we draw attention to an al- 
ternative interpretation in which no planetesimal belt is 
required, but rather some of the disk structure is associ- 



ated with grain destruction or alteration by sublimation 
or annealing. 

4. By comparing the imaging data with that acquired 
in 2003, we confirmed the existence of the dusty clump 
at 52 AU in the SW disk. We speculate that the clumps 
projected spatial displacement of ^2.0 AU, a 3.6-er result, 
between two epochs separated by seven years may be due 
to the Keplerian motion of the clump at an orbital radius 
of 54.3±f;§ AU. 
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